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ABSTRACT 



The persistent increases in spin-down rate {offsets) seen to accompany 
glitches in the Crab and other pulsars suggest increases in the spin-down 
torque. We interpret these offsets as due to starquakes occurring as the star 
spins down and the rigid crust becomes less oblate. We study the evolution 
of strain in the crust, the initiation of starquakes, and possible consequences 
for magnetic field and torque evolution. Crust cracking occurs as equatorial 
material shears under the compressive forces arising from the star's decreasing 
circumference, and matter moves to higher latitudes along a fault inclined to 
the equator. A starquake is most likely to originate near one of the two points 
on the rotational equator farthest from the magnetic poles. The material breaks 
along a fault approximately aligned with the magnetic poles. We suggest that 
the observed offsets come about when a starquake perturbs the star's mass 
distribution, producing a misalignment of the angular momentum and spin 
axes. Subsequently, damped precession to a new rotational state increases 
the angle a between the rotation and magnetic axes. The resulting increase 
in external torque appears as a permanent increase in the spin-down rate. 
Repeated starquakes would continue to increase a, making the pulsar more of 
an orthogonal rotator. 



Subject headings: dense matter — magnetic fields — stars: magnetic fields — 
stars: neutron — pulsars: individual (Crab) 
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1. Introduction 

The magnetic braking torque acting on an isolated neutron star would be steady in the 
absence of abrupt changes to the star's magnetic configuration. Most pulsars, however, do 
not slow in a regular fashion, but suffer variations in their spin rates in the form of glitches 
and timing noise. Perhaps the most striking aspect of spin evolution is the persistent 
increases in spin-down rate that accompany glitches in the Crab pulsar ([Lyne, Pritchard 



Smith 19951 ), PSR 0355+54 ( |Lyne 1987| ) and PSR 1830-08 QShemar fc Lyne 1996| ); see 



Table 1. In the Crab, these permanent offsets involve fractional changes in the spin down 
rate of ~ 10^^ - 10^^ (see Fig. 1). PSR 1830-08 has exhibited one offset of 8 x 10~\ and a 
persistent offset might have followed the large 1986 glitch of PSR 0355+54. Similar offsets 
might also occur in the Vela pulsar ( [Link fc Epstein 1997|) , giving rise to the small braking 



index of 1.4 ± 0.2 reported by Lyne et al. (1996). 

It is striking that all observed offsets are of the same sign, and correspond to 
increases in the spin-down rate. One interpretation of this phenomenon is that glitches are 
accompanied by sudden and permanent increases in the external torque acting upon the 
star (pijink, Epstein fc Baym 1992| ; P^ink fc Epstein 1997]) . Such torque changes could occur 



if either the direction or magnitude of the star's magnetic moment changes. Starquakes, 
occurring as the star spins down ( |Baym fc Pines 1971 ; Ruderman 1976 ) would affect the 



external torque if they change the orientation of the magnetic moment with respect to 
the spin axis. If structural relaxation occurs asymmetrically about the rotation axis, due 
perhaps to magnetic stresses or asymmetric material properties, the star's spin and angular 
momentum vectors would become misaligned. As the star precesses and relaxes to a new 
rotational state, the magnetic moment would assume a new orientation with respect to 
the rotation axis, leading to a change in the external torque. In this paper, we study how 
the rigid neutron star crust relaxes its structure, and consider possible consequences for 
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evolution of the magnetic field and torque. 



Alpar and Pines (1993; also [Alpar et al. 1996]) have suggested that the Crab's offsets 



result from a reduction in the moment of inertia on which the external torque acts. Such 
a change could occur by either a structural change of the star, e.g., the star becomes 
less oblate, or through a decoupling of a portion of the star's liquid interior from the 
external torque. If the moment of inertia decreases through structural readjustment, 
to conserve angular momentum, the star would always spin more rapidly than had the 
glitch not occurred ([Link, Epstein fc Baym 1992| ; [Link fc Epstein 1997| ). The large offsets 
following the 1975 and 1989 glitches in the Crab, however, eventually caused the star to 
spin less rapidly than had the glitch not occurred ( |Lyne, Pritchard, fc Smith 1993 ); the 



observed offsets, therefore, cannot be due solely to structural readjustments. In principle, 
the Crab's spin-down rate offsets could be due to decoupling of a portion of the star's 
superfluid interior from the external torque ( [Alpar fc Pines 1993| ; [Alpar et al. 1996]) , though 



quantitative agreement of this model with the data has yet to be demonstrated. Torque 
increases associated with the surface field structure appears to be the most straightforward 
explanation. 

A neutron star relaxes its oblateness as it spins down, moving equatorial material 
toward the rotation axis and reducing the equatorial circumference. If the stellar crust is 
brittle, stresses lead to starquakes as the yield strength of the crustal material is exceeded 
( [Baym fc Pines 1971[ ).f[ The actual response of solid neutron-star matter at high pressure 



to shear is not well-understood. Here we assume that the neutron star crust is brittle and 
explore the consequences of this assumption.^ 



^Superfluid stresses in the crust ( [Ruder man 1976| ) could also drive crust cracking, but 
will not be considered here. 

^ Known materials exhibit ductile rather than brittle behavior under pressures comparable 
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2. Growth of Strain in a Spinning-down Star 

To study the quahtative features of the growth of strain in a neutron star as it 
spins-down under its external torque, we model the star as an incompressible, homogeneous, 
self-gravitating elastic sphere of constant shear modulus. The non-rotating configuration 
is a sphere of radius R. We take the initial configuration to be unstrained with rotation 
frequency Q, and investigate how the star deforms as the spin rate changes by —6Q 
where < ^ The initial rotation produces an equatorial bulge of relative height 
~ {RVt/vkY ^ 1) where Vk = {2GM / R)^^"^ is the Keplerian velocity near the stellar surface. 
A change —5VL in the rotation rate produces a small decrease in the equatorial bulge of 
~ VtdVLR^/vf.. As the star deforms, matter originally at r is displaced to r + 'u(r), where 
u{r) is the displacement field. The displacement field can be found using elasticity theory 
(see, e.g.^ [Love 192"0| ; [Landau fc Lifshitz 1959|) . In a spherical coordinate system centered 
on the star, with 6* = at the star's north pole, the displacement field is (see, also, [Baym fc 
Pines 1971[) 



u^{^r,e) = ^M^Wr^-^) (l-3cos^^ 



VtdVtR^ 



U0(r,6) = K — r (5r — 8) sin 6* cos 6*. 

4 ^ ' 



Here r and u are in units of the stellar radius, and we have neglected terms of order 
{ct/vkY'^ 10"^ where q is the transverse sound speed. The displacement field is illustrated 
in Fig. 2. 



to their material shear moduli ( puba et al. 1990[) . Nevertheless, deep-focus earthquakes are 



known to originate from regions of very high pressure ( [Green fc Houston 1995[ ) . These faults 
are thought to be facilitated by densification phase changes; small regions of higher density 
nucleate as the material is stressed, and act as a lubricant for shearing motion. Analogous 
processes might occur in the high-pressure material of the neutron star. 
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The local distortion of the solid is described by the strain tensor 

In a local coordinate system in which this matrix is diagonal, the eigenvalues represent 
compression (e^ < 0) or dilation (e, > 0) along the respective axes. Suppose the largest 
and smallest eigenvalues are ei and e^, respectively. The strain angle is e/ — e^, and the 
plane containing the corresponding principal axes is the stress plane. For the uniform 
sphere, we find that the stress plane is tangential to the stellar surface near the equator 
with compression in the azimuthal direction and dilation toward the rotational poles (see 
Fig. 3). Near the poles, the stress plane switches to the radial-azimuthal plane, i.e., along 
cones. These results are consistent with those of Baym and Pines (1971), who calculated 
the strain field as a sphere at zero pressure spins down. 



3. Star quakes 

At zero pressure, a stressed element of matter can break due to shear or tension. At the 
pressure of the neutron star crust {p ^ pCj), however, only the former process can occur. 
Compressed material strained at a critical angle shears along a plane as shown in Fig. 4 
( preen fc Houston 1995|) . The angle this plane takes with respect to the plane of maximum 



compression depends on the material's internal friction, and is in the range 30 — 45°. The 
critical strain angle of neutron star material is not accurately known. Highly-compressed 
terrestrial minerals can support strain angles of >0.01 before breaking (pillings 1972| ; Duba 
et al. 19901) . 



The geometry of a starquake is independent of the value of the critical strain angle 
for neutron star matter. Strain builds as the star slows and reduces its circumference. For 
the homogeneous sphere, the maximum strain is reached on the equator, as each material 
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element is compressed by its neighbors on the equator (see Fig. 3). When the strain 
angle for an element reaches a critical value, the material breaks along a fault that takes 
an angle of 30 — 45° with respect to the equator and propagates to higher latitudes (see 
Fig. 5). Movement of material from the equatorial bulge to higher latitudes decreases the 
star's oblateness. At higher latitudes the strain angle diminishes and the fault terminates 
( [Hertzberg 1996| ). This picture of starquakes differs fundamentally from that of Baym & 
Pines (1971) who considered the distortion in a sphere at zero pressure. They found that 
the material breaks as tensile forces pull apart material in the equatorial plane. In our 
analysis, however, the enormous pressure in neutron star matter prevents tensile rupturing. 

For an axially-symmetric star, the crust is equally likely to begin breaking at every 
point on the equator. The rotation-induced displacements of the star's crust distort the 
initial magnetic field B{r) that is anchored in the stellar crust thus generating magnetic 
stresses that affect the development of strain. The strength of the magnetic stress relative 
to the shear stress is ~ /3 = {vA/ct)'^, where va is the Alfven velocity. Through most of the 
crust the induced magnetic stresses are significantly smaller than the material shear stresses 
and can be treated perturbatively. For example, at moderate densities in the outer crust 
we have (3 = 8 x 10-^(S/10^2q)2(^/^qio g ^^-a^-i^^^/^O^cm s~^y^, while the magnetic 
stresses can exceed the material ones at low densities, below ~ lO^gcm"^. 

The strain field that develops in magnetized material could be obtained by solving 
the elasticity equations and Maxwell's equations self-consistently, subject to boundary 
conditions at the star's surface. To illustrate how the star's magnetic field determines where 
starquakes originate, we consider the simpler problem of how the presence of a magnetic 
field affects the material strain in an infinite medium. In the absence of the field, the 
equilibrium state is given by 
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where Fj represents gravitational and centrifugal forces. In the presence of a magnetic field, 
the stress tensor changes by Suij. The displacement of the material induces a change in 
the field with a corresponding change in the Maxwell stress tensor T^j of 5Tij. We take 

— * 

the unperturbed field B to be dipolar, which gives no force on the undistorted matter, ie., 
dTij/dxi — 0. The new equilibrium is given by 

±.{Sa,j + ST,j) ^ 0. (4) 

For an infinite medium (one in which boundary effects are unimportant), the solution is 
Saij = —STij. The magnetic field induces a change 5u in the displacement field. For a 
field that is sufficiently small that this correction to the displacement field is small, the 
perturbation to the Maxwell stress tensor is 

5Tij = [SBiBj + BiSBj - B^SB^Sij] , (5) 

with 

5B^Vx(uxB). (6) 

Here u is the displacement field in the absence of the field. For incompressible matter, the 
correction 5uij to the strain tensor is 

where /i is the shear modulus. The eigenvalues of this tensor give the correction to the 
strain angle due to the induced magnetic forces. Using the displacement field of eq. [1] in 
eq. [6], we obtain an estimate of the corrections to the strain angle on the star's equator; the 
result is shown in Fig. 6. Near the magnetic poles, the field gives the material extra rigidity, 
and reduces the strain angle. Consequently, as the star spins down, the strain increases 
the fastest at the two equatorial points farthest from the poles. These results indicate that 
a starquake will begin at one of these two points when critical strain is reached. A more 
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complete treatment would account for magnetic stresses at the stellar surface; we do not 
expect treatment of the boundary to alter our qualitative conclusions. 

Isotropic material under compressive stress is as likely to break along the plane 
indicated in Fig. 4 as along its complementary plane. In a neutron star, however, the 
magnetic field gives the material extra rigidity along shear planes that cross field lines. 
Cracking is thus favored along the plane that is most parallel to the field (fault / in Fig. 5). 

4. Effects on the Spin-down Torque 

In the picture of starquakes described above, a slowing star reduces its equatorial 
circumference by shearing material across the equator, moving the material along the fault 
to higher latitudes (see Fig. 5). We now discuss how this crustal motion could affect 
the spin-down torque acting on the star by producing a misalignment of the star's spin 
and angular momentum vectors. The star then undergoes damped precession to a new 
rotational state.0 

Quake-induced mass motion produces a small change in the orientation of the principal 
axes of the star's moment of inertia. The material motions create "mountains" (~ 100 
/im high) at higher latitudes, breaking the axial symmetry of the star's mass distribution. 
Suppose that before the starquake the largest principal moment of inertia corresponds to 
a principal axis aligned with the rotation axis. As the mass motion occurs, the principal 
axis of inertia shifts away from the magnetic axis by Aa, pointing in a new direction ei 
(fixed in the star, see Fig. 5). The star's rotation and angular momentum vectors are now 

^Crust motion would also change the external torque by moving the magnetic poles with 
respect to the rotation axis. The associated change in a, however, is only ~ AR/R ~ 10^^, 
too small to account for the largest persistent shifts seen in the Crab. 
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misaligned and the rotation axis precesses about the angular momentum axis ( ^hahani 



1977|) . Dissipative processes, such as the interaction of superfluid vortices with the nuclei 



of the inner crust, damp the precessional motion (|Sedrakian, Wasserman fc Cordes 1998|) , 



bringing the star to a new stable equilibrium in which the principal axis of inertia is again 
aligned with the star's angular momentum. As the precessions damps, the angle between 
the star's rotation and magnetic axes a increases by Aa. In some models of pulsar spin 
down, an increase in a leads to an increase in the spin-down torque. In the magnetic dipole 
model for pulsar slowdown, for example, an increase by Aa gives a permanent increase in 
the spin-down rate of Afi/fi = 2Aa/ tana. We suggest that the Crab's offsets arise from 
starquake-induced changes of order Aa ~ 10""^ ( [Link fc Epstein 1997|) . 



The following dimensional argument gives an estimate of the magnitude of the shift of 
the principal axes. Suppose that between starquakes the star spins down by 5VL and reduces 
its equatorial circumference by Ai?. As the starquake occurs, crust material suddenly 
moves along the fault by a comparable amount. From eq. |l|], we estimate 



Ai? ^ i^T^R^nsn _ ^ f Q y / V, Y'( ]( W y 

R 3vl V200 rads-V VlO^o cm s-V \UyiJ\10^jrJ '^^ 

where tg is the time interval between large quakes and tage is the spin-down age. Assuming 
that quakes are coincident with glitches, tg ~ 14 yr for the Crab pulsar. Matter accumulates 
in mountains over a length scale of lateral extent comparable to the fault length L. We 
expect L to be roughly the length scale characterizing the strain field, of order R. Upon 
formation of mountains, ei changes its angle with respect to the rotation axis by an amount 
of order the ratio of the off-diagonal elements of the moment of inertia tensor to the bulge 
moment of inertia; the shift is of order Aa ~ {M^/ Mb){L/ R), where is the mass of a 
mountain and Mb is the mass in the equatorial bulge. The mass in the bulge is of order 
Mb ~ [RQ/vkYMc, where Mc is the mass of the crust. The fraction of crust mass that 
accumulates in a mountain is roughly the fraction of the star's area that moves with the 
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fault, i.e., Mm ~ {ARL/R'^)Mc. Hence, we obtain an upper limit of 

3 R^Q \rJ VMyry VlO^yry ^' 

For the Crab pulsar, Rankin (1990) obtains a = 86° from radio data, and Yadigaroglu 
& Romani (1995) find a = 80° from their analysis of gamma-ray and radio data. The 
corresponding fractional change in spin-down rate (=2Aa/ tana) is ^10^^ for a = 80°, 
large enough to explain the observed persistent shifts; a = 86° gives a limit of 2 x 10^^. 
We cannot yet apply eq. |^ to PSRs 0355+54 and 1830-08; only one large glitch has been 
observed in each of these pulsars, so the quake intervals and typical offset magnitudes are 
unknown. Moreover, the magnitude of the permanent offset in PSR 0355+54, if any, is 
unknown since the spin-down rate has continued to changed following the large glitch of 
1986. 



5. Summary 

The persistent increases in spin-down rate {offsets) seen to accompany glitches in the 
Crab and other pulsars suggest sudden increases in the spin-down torque. Starquakes 
occurring as a neutron star spins-down and readjusts its structure affect the spin-down 
torque exerted on the star by changing the magnetic field geometry or orientation. In this 
paper we have examined the evolution of strain in the crust of a spinning-down neutron 
star and the initiation of starquakes as the material reaches critical strain. Crust cracking 
occurs as equatorial material shears under the compressive forces arising from the star's 
decreasing circumference. The star decreases its oblateness as matter is moved to higher 
latitudes along a fault inclined to the equator. Magnetic stresses suppress shearing near the 
magnetic poles, especially shearing motions across the field lines. Starquakes are thus most 
likely to originate near the two points on the equator farthest from the magnetic poles and 
propagate toward the magnetic poles. 
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Starquake-induced misalignment of the star's angular momentum and spin, associated 
with glitches, is a possible explanation for the spin-down offsets seen in the Crab pulsar. 
Following the misalignment, damped precession to a state of larger angle a between the 
magnetic and rotation axes could increase the external torque, giving a permanent increase 
in the spin-down rate. Repeated starquakes would continue to increase a, making the 
pulsar more of an orthogonal rotator. 

We thank A. Ohnto, P. Richards and D. Longcope for valuable discussions. This work 
was performed under the auspices of the U.S. Department of Energy, and was supported in 
part by NASA EPSCoR Grant #291748, and by IGPP at LANL. 
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TABLE 1 



PSR 


age (yr) 


glitch year 


AQ/Q 




permanent offset AQ/ fi| 


Crab^ 


10^ 


1969 


4 X 10" 


-9 


-4 X 10-^ 






1975 


4 X 10" 


-8 


-2 X 10"^ 






1986 


4 X 10" 


-9 


-2 X 10^^ 






1989 


9 X 10" 


-8 


-4 X 10-^ 


0355+542 


6 X 10^ 


1986^ 


4 X 10- 


-6 


-3 X 10-3? 


1830-082 


2 X 10^ 


1990 


2 X 10" 


-6 


-8 X 10"^ 



^Lyne, Pritcfiard & Smitfi 1993; ^Sfiemar & Lyne 1996 

^Tlie magnitude of this offset is quite uncertain, as the pulsar has not had time to completely 
recover from the glitch. 
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Fig. 1. — Twenty- five years of spin history of the Crab pulsar showing apparent torque 
increases (adapted from Lyne, Pritchard, & Smith 1993). Shown are spin rate residuals 
relative to a model for data prior to the first glitch. Glitches are indicated by arrows. 
Following each glitch, the pulsar acquired a greater spin-down rate than before the glitch. 
The small glitch of 1986 was also followed by a small offset, shown by the dashed fine. These 
offsets appear to be permanent and cumulative. 

Fig. 2. — The displacement field in a spinning-down neutron star. A cross-section through 
the center of the star is shown. The equatorial diameter decreases, while the polar diameter 
increases. 

Fig. 3. — Strain eigenvalues and strain angle vs. latitude 9 at the stellar surface for B — 
(arbitrary units). The short-dashed curves give the strain angle; the other curves are the 
eigenvalues. 

Fig. 4. — Breaking of a compressed block of matter. A block under horizontal compression 
and vertical tension shears along a plane as shown when the critical strain is reached. 
Shearing along a complementary plane fiipped over with respect to the plane shown is equally 
likely for isotropic material. The y — z plane is the stress plane. 

Fig. 5. — A starquake. In the absence of a magnetic field, the material is equally hkely to 
begin breaking along faults / and /', anywhere on the equator. In the presence of magnetic 
stresses, fault / is more likely, creating "mountains" (indicated by the snow-capped peaks) 
and shifting the largest principal axis of inertia to a new direction ei (fixed in the star). 

Fig. 6. — The strain angle of surface equatorial material in the presence of a magnetic field 
(arbitrary units). The curves correspond to different values of the angle a between the 
magnetic and rotation axes. The magnetic poles are at azimuthal angles 7r/2 and 37r/2. 
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